General Solution of EM Wave Propagation in Anisotropic Media 
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When anisotropy is involved, the wave equation becomes simultaneous partial differential equa- 
tions that are not easily solved. Moreover, when the anisotropy occurs due to both permittivity and 
permeability, these equations are insolvable without a numerical or an approximate method. The 
problem is essentially due to the fact neither e nor jj, can be extracted from the curl term, when they 
are in it. The terms V x E (or H) and V x eE (or /iH) are practically independent variables, and 
E and H are coupled to each other. However, if Maxwell's equations are manipulated in a different 
way, new wave equations are obtained. The obtained equations can be applied in anisotropic, as 
well as isotropic, cases. In addition, E and H are decoupled in the new equations, so the equations 
can be solved analytically by using tensor Green's functions. 
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I. INTRODUCTION 

As the importance of anisotropic devices has increased 
in many fields of optics and microwaves, wave propaga- 
tion in anisotropic media has been widely studied over 
the last decades The anisotropic nature basically 
stems from the polarization or magnetization that 
can occur in materials when external fields pass by. 
Generally, certain axis components of the E and the H 
fields are influenced by other axis components and by 
those of the same axis. This is why matrices are involved 
in e and \i. Therefore, Maxwell's equations and the wave 
equations are also represented in matrix form. Mathe- 
matically the electric field E and the magnetic field H are 
not only vectors, but also rank-1 tensors, which implies 
that they obey the set of rules of coordinates transforma- 
tion. This also implies that e and p, are tensors of rank-2. 

In ordinary homogeneous isotropic media, a wave 
equation with a source term is solved with an ordinary 
adaptation of Green's functions. However, in anisotropic 
media, the equations become linear simultaneous partial 
differential equations (PDE). These equations contain 
all information pertaining to the anisotropic prop- 
erties. As they are linear, Green's functions can be 
considered easily with the source terms. These types 
of solutions with Green's functions have been studied 
under the name of "dyadic Green's functions" over the 
last few decades. Most of published papers presented 
the special properties of applications or methods of ob- 
taining Green's functions in a specified set of coordinates. 



This paper presents a general methodology of solving 
the wave propagation equation in an anisotropic envi- 
ronment to obtain the E and the H fields. The tensor 
Green's functions GEij and GHij for an electric and a 
magnetic field are used, as these equations are linear. 
Generally for an unbounded case, a Fourier-transformed 
Green's function is useful. It changes differential equa- 
tions into algebraic equations in k-space. Therefore, 
finding the Fourier-transformed Green's functions is the 
core procedure in solving the wave equations. 



In Section II, the anisotropic characteristics are only 
concerned with the permittivity e, which becomes a ma- 
trix. The corresponding Fourier-transformed Green's 
functions for this case are easily obtained. Therefore the 
equation system is solvable. The Section III addresses a 
case in which anisotropy is described only with the matrix 
permeability \x. The mathematical procedure in this case 
is similar to that discussed in Section II. The subsequent 
section, Section IV, is an arbitrary case. The anisotropy 
comes from both the permittivity and the permeability. 
The equations appear to be insolvable in an analytical 
sense. Other algebraic manipulations of Maxwell's equa- 
tions lead to new wave equations, in which E and H arc 
decoupled. Analytic solutions are also possible with the 
use of Green's functions. The last section contains the 
discussion and conclusion. 



II. ANISOTROPY FROM PERMITTIVITY 



'Electronic address: arameas@kaist.ac.kr 
t Electronic address: stlee@semyung.ac.kr 



The typical method for deriving of the wave equation 
starts from 
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V x V x E = -g-(V x juH), 



(1) 



VxVxH = VxJ + —(V x eE). (2) 

at 

In an ordinary homogeneous isotropic case, the constants 
e and p are extracted from the curl term, and the equa- 
tions lead to wave equations. In this section, it is as- 
sumed that anisotropy exists due to the permittivity. Ad- 
ditionally, e is a matrix, but p is a scalar number. Thus, 
Eq.(l) is considered as an isotropic case. On the other 
hand, Eq.(2) has a problem in the curl term on the right 
side. The matrix e cannot come out of the curl term, and 
Maxwell's 3rd equation cannot be applied. Therefore, the 
equations become 



V x V x E - 



1 <9 2 E _ 
d 



V^H = -V X J - eo^(V x e r E). 



(3) 
(4) 



Here the first term of Eq.(3) was not expanded to 
V(V • E) — V 2 E, as the divergence of E cannot be 
written as e~ 1 p if using Maxwell's equation. When e is a 
matrix, the correct form of the first Maxwell's equation 
is V • eE = p, instead of V • E = (T 1 p. The e cannot 
move out from the inside of divergence. 



Except for the fact that e is matrix, Eq.(3) resembles 
an isotropic case. It is and equation for E, but it is 
important to note that the e in the curl of Eq.(4) cannot 
come to the front of the curl operator, as a constant 
number does. If the matrix satisfying 4-VxE = Vxe-E 
is found, Maxwell's equation becomes applicable. How- 
ever there is unfortunately no matrix like A. Therefore, 
eE in the curl term serves as an independent variable 
in a practical sense. Thus, E and H are coupled in the 
equation, which is not solvable by itself. Nevertheless 
Eq.(3) can be analytically solved using tensor Green's 
functions. The tensor Green's function is explained in 
the appendix to check if it satisfies the Green's function 
condition and whether or not it solves the problem. 

After obtaining the solution E, the result is substituted 
into Eq.(4). A Solution for H then becomes possible 
with the Green's function technique as well. In this case, 
the term (V x e r E) on the right side of Eq.(4) becomes 
a part of the source terms, as E has already been solved 
from Eq.(3). 

To solve the equations, it is convenient to assume that 
the time dependency is harmonic, or e lLJt , and to use the 
following notations for the source terms of Eqs.(3) and 
(4): 

U(r) = -iujp p r J, (5) 
V(r) = -V x J - iwe V x (e r E). (6) 



The corresponding Green's functions satisfies the follow- 
ing conditions: 



V x V x GEW(r,r') - 

= <5(r - r'), 
V 2 Gff (1) (r,r') = 6{r- 



k 2 GE^(i 



(7) 
(8) 



GE and GH are the Green's functions for an electric and 
a magnetic field, respectively. The superscript (1) in GE 
and GH serves to distinguish the anisotropy that occurs. 
1 denotes this for e, 2 for p, and 3 denotes this for both e 
and pL. The Fourier transform is effective for an unbound 
case: 



G£ (1) (r,r') = /// gE {1) (k) e^ k (r - r Vk. (9) 



Inserting the Fourier-transformed function in Eq.(9) 
into Eq.(7), we obtain the algebraic equations 
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£«(k) 



1 



((/c 2 I-k®k)-fc 2 ^ r e r )' 



(10) 



Here, <S) refers to the direct product, and I is an identity 
matrix. The above equation for gE^(k.) is the main 
step involved in the solution. Cottis et al. calculated 
this Green's functions in cylindrical coordinates [2j. The 
GE for the Green's functions could be calculated by 
inserting Eq.(10) into Eq.(9). 

The form of Eq.(8) is identical to the Poisson equa- 
tion in electrostatics. The following equation is, thus, 
applicable: 



47r r - 



(11) 



The solutions can then be written as follows according to 
the usual Green's functions method: 



E(r) = JJJ GE (1) (r,r')U(r')dV, ( 12 ) 
H(r) = [[[ GH^(r,r')V(r')d 3 r'. (13) 



Looking at Eq.(3), (— V 2 E) is concealed in (Vx VxE), 
and the equation is a type of wave equation. As e r is a 
matrix, E\, E% and £3 are coupled to each other. There- 
fore, this equation is fundamentally simultaneous PDE. If 
e r is diagonalized, Ei 's becomes decoupled after replac- 
ing (V x V x E) by V(V • E) - V 2 E . If the eigenvalues 
are Ai,A2 and A3, the wave vector changes from fej to 
ki\f\~i on the principal axis. Therefore, we can compute 
its refractive index or change of the wave velocity along 
the corresponding axis. 
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III. ANISOTROPY FROM PERMEABILITY 

When anisotropy comes only from \i r , fir becomes a 
matrix, and e r remains a scalar number. The basic equa- 
tions are as follows: 

V x V x E = -iw/xo(V x (/i r H)), (14) 
V 2 H + fc 2 e r /x r H = -V x J. (15) 

As in the previous section, V x (/i r H) in Eq.(14) makes 
the problem complex. This equation cannot be solved di- 
rectly. However, Eq.(15) is merely Helmholtz-type equa- 
tions with the source term, — V x J. The calculation is 
carried out by using a Fourier transform, as in the pre- 
ceding section. The Green's functions for Eqs.(14) and 
(15) are as follows: 

VxVxG£ (2) (r,r')=i(r-r'), (16) 
V 2 Gff (2) (r,r') + fc 2 e rAV G.ff (2) (r,r') 

= <5(r-r'). (17) 

The introduction of the Green's functions GE and GH 
proceeds identically as it did before. Eq.(17) is a general 
Helmholtz-type equation of the type studied by many 
authors [3j. The Fourier-transformed Green's function 
gE^ 2 \k) and gH^ 2 \k) are defined as follows: 

GE™ ( r , r ') = JJJ gE W (k) e - lk < r - r Vk , (18) 

Gff( 2 »(r,r') = J I J .gtf( 2 )(k)e-* k -( r - r Vk . (19) 

Hence, the algebraic form of gE^ 2 \k), gH^ 2 \k) is ob- 
tained through insertions into Eq.(16) and (17): 



(fc 2 I-k®k) ' 



(20) 
(21) 



where the denominator is also a matrix showing a tensor 
property. At this point, the solution for H is possible 
using GH( 2 \r, r'): 



H(r) 



GiJ (2) (r,r') V(r')dV. 



(22) 



The source term in the above equation is V(r) = —V x 
J. The next step is to find the solution. As the function 
for H, has already been solved, the source term of Eq.(14) 
becomes U(r) = — zw/io(V x (/i r H)). The solution E is 
then 



E(r) = /// G£ (2) (r,r') U(r')dV. 



Now, the fields E and H are solved when the 
anisotropic characteristic comes from either e r or fir- 
Problems were noted in dealing with the curl term in- 
cluding an anisotropic factor. However, this difficulty is 
circumvented by solving the other equation first and by 
placing the result into the curl term, which causes prob- 
lem. The solutions can be obtained, but the equations 
loose their original wave shapes. This does not mean 
they are not a wave, but the analytical wave nature is 
not directly observable. 



IV. ANISOTROPY FROM PERMITTIVITY 
AND PERMEABILITY 



When anisotropy occurs due to both e and fi, the wave 
equations become more complicated to solve. This diffi- 
culty arises in the same way. The terms e r • E and fi r ■ H 
are merely a linear combination of the original E and H, 
but when contained inside the curl term ( V x ) , they act 
as independent variables. For example, when there is a 
term (V x E), we can replace it by (— ^) using Maxwell's 
equations. However, this replacement cannot be applied 
to the term (V x e r E) . It is practically another unknown 
to extent that E is solved. The basic wave equations take 
the following forms: 



Vx(VxE) = -iw/ioV x (/i r H), (24) 
V 2 H = -V x J - iwe V x (e r E). (25) 

The goal here is clearly to find E and H, but there are 
two more terms that invoke problems, (V x fi r H) and 
(V x e r E). Therefore, it is impossible to obtain the ana- 
lytical solution simultaneously from the wave equations. 
A better approach is to go back to the original Maxwell's 
equations and derive new equations instead of relying on 
the original wave equations: 



VxE = —pa 



d(/irH) 



V x H = J + e 



at ' 

d(e r F,) 

at ' 



(26) 
(27) 



The magnetic field H is written from Eq.(26) as 

H = / — Hr 1 ■ (V x E) dt. (28) 
J Mo 

Inserting the above equations into Eq.(27) and differen- 
tiating with respect to time gives 



V x (h; 1 ■ (V x E)) + 1 



-Mo 



dJ 

at 



(29) 



(23) 



This equation restores the original wave equations when 
the media is isotropic. Moreover, it is an equation for E 
that is decoupled from H. In fact, this equation holds in 
any case, regardless of the existence of anisotropy. If the 
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time dependency in E and J are assumed to be harmonic 
or e Mt , the above equations read. 

V x (/i" 1 • (V x E)) - kl e,.E = -icj/i J (30) 

As the equation is linear, the tensor Green's function 

Vx(/*- 1 -(VxG£( 3 >(r,r')))-fc§ e r G£ (3) (r, r') = J(r-r') 

(31) 

is also possible. Assuming that gE^ (r, r') is Fourier- 
transformed in the same way to the Eq.(18) in the last 
section, the spatial Green's function 

G£ (3) (r,r') = [[[ <?£ (3) (k)e- ,;k - (r ~ r Vk. (32) 



Eq.(31) provides the algebraic relationship for gE^: 

-kx/i^xfkx gE^) - kle r gE^) = 1. (33) 

By applying the following matrix k representing the wave 
vector k Hi, 



(34) 






-k z 


Ky 


h 





-k 


— ky 








Eq.(33) above becomes 
<?£ (3) (k) = 



(35) 



Here, k is defined in rectangular coordinates. This 
implies that k is coordinates dependent. Now, <?£^ 3 )(k) 
is obtained, and the calculation of GE^ and E(r) are 
straightforward. 



G£ (3) (r,r') = JJJ <?£ (3) (k) e- jk ( r - r Vk , (36) 
E(r) = /// G£ (3) (r,r') U(r')dV , (37) 



where U(r) is the right term of Eq.(30) as a source. The 
solution of H(r) is obtained directly from Eq.(28) by in- 
serting E. In this section, a classical wave equation itself 
was not considered. Instead, a derivation of new solvable 
wave equations from Maxwell's equations was presented. 



V. DISCUSSION AND CONCLUSION 

When anisotropy occurs due to either permittivity 
or permeability, the wave equation is solved by using 
Green's functions. However, in an arbitrary case in 
which anisotropy occurs in both e r and n r , the equation 
is not substantially solvable in an analytic sense. The 
usual means is to resort to a numerical calculation. 



However, there is another method that gradually 
approaches solutions by iteration. This does not need to 
be a numerical calculation. It can be performed analyt- 
ically. Either e r or u r becomes a scalar number; then, 
the whole system can be solved as before. Therefore, the 
solution of the equations becomes feasible by replacing 
e r with a scalar number. The solution obtained in this 
way is the first trial solution. The equations are then 
solved again by inserting the first solutions into the curl 
term, but from this time, the original matrix value for 
e r is used instead of the first number that was chosen 
at the beginning. An improved second solution for E 
and H is then calculated. If this iteration continues, 
the results approach some converging functions, which 
are supposed to be solutions to the equations. This 
procedure can be carried out until difference between 
the (n — l) th solution and the n th solution becomes 
less than a prescribed level. This process can be done by 
using a symbolic calculation [5j. However, it is merely 
an approximate method, regardless of whether it is a 
numeric or a symbolic calculation. 

Instead of solving the original equations, we derived 
new wave equations that were especially useful for 
an anisotropic problem via a slight manipulation of 
Maxwell's equations. The new equation is identical to 
the original wave equation when there is no anisotropic 
characteristic. The fundamental advantage of this new 
equation is the fact that E and H are decoupled, unlike 
in the original equations, making it possible to find 
analytical propagators. 

In an actual calculation, there is one factor to consider; 
i.e. Eq.(35) becomes singular in the isotropic case. The 
determinant of the denominator becomes zero in this 
special case when e r — fj. r — scalar number. In such 
an isotropic case, the solution can be obtained in an 
ordinary manner. 

Eq.(29) was derived for an electric field. An equation 
for the magnetic field can also be obtained in a similar 
way. The result is as follows: 

Vx^VxHJ + ifr^Vx^J). (38) 

Solving this equation is identical to the procedure for E 
given in the last section. 

The new wave equations can be justified in several 
ways. The best approach is to compare the experimen- 
tal values with the calculated result for the new equa- 
tions. However, many aspects can be examined analyti- 
cally by comparing the results of the two systems of the 
equations. The first comparison is that the two equa- 
tion sets are identical in an isotropic case, these can 
be seen immediately by inspection, because there is no 
matrix. For anisotropic cases, the Fourier-transformed 
Green's functions can be checked whether they are iden- 
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tical or not. Cottis, Vazouras and Spyrou calculated the 
Fourier-transformed dyadic Green's functions in a dielec- 
tric anisotropy, and their result is the same as Eq.(lO). 
The result of the new equation for an identical case is 
given below according to Eq.(35): 



jE^(k) 



-1 



(39) 



Given that there is no anisotropy in /li, the identity 
matrix is substituted in place of \i~ x . The same re- 
sult is obtained using the relationship k ■ k = k<x>k— k 2 I. 

The other comparison is for the case of magnetic 
anisotropy. In this case e r = I. The result of the Fourier- 
transformed Green's functions of the original equation is 
given in Eq.(21). It can be compared with the result of 
Eq.(38), which gives 



gH 



(3) 



-1 



k €>p 1 k ~\~ /Cq fjjy 

k I k ~\~ f-if 



k ® k - k 2 I + kl n r 



(40) 



However, the k ® k term becomes zero when it is 
applied to a magnetic field. Hence, the result of the 
original equation is identical to that of the new equations. 

For more verification, the plane wave E = 
E exp[i(k ■ r — cot)] can be considered in an anisotropic 
medium. When the relative permittivity tensor is given 
in terms of refractive indices as 




$ 



ni 



(41) 



the refractive indices then satisfy the following well- 
known condition: 



-n 2 (k 2 k 



Th k^£ ky 

-n 2 (k 2 x + k 2 ) 



*k-rrk?_ 



n z — n (k x k y ) 



(42) 

where the conventional notations n = £;k and fcj = p^j- 
are used. Exactly the same result is also obtained from 
the new wave equation, Eq.(29). This is actually a 
natural consequences, because they are based on the 
same mathematical ground. 

As for advantages of the new equations, it is clear that 
original equations do not allow an analytical approach 
in the general case whereas the new equations give 
analytical propagators (Fourier-transformed Green's 
functions). The point is whether or not they are 



decoupled. This is a considerable difference between the 
two systems, although the new equations need numerical 
calculations, for example integration to get GE(r,r'), 
GH(r, r') or the final answers E and H. The existence 
of propagators is supposed to bring a non-negligible 
amount of code-saving effect. 

As a result, an equation set Eqs.(29) and (38), is ob- 
tained and describes wave propagation in anisotropic me- 
dia. This shows that the two systems are commutable 
and equivalent. The new wave equation is expected to 
be useful in analytical research of anisotropic properties 
beyond what a numerical approach can address. 



Appendix A: Proof of the Tensor Green's function 

As Eq.(3) is linear, it is natural to consider Green's 
functions. However, the equation is a multi-linear simul- 
taneous PDE, and can verify if the use of Green's func- 
tions is proper. If Green's functions work, the following 
type of solutions can be assumed as valid: 



E,; 



GEijVjdr' 



(Al) 



The functions GEij are devised to play the roles of in- 
fluence functions to generate the fields due to p and J. At 
this stage, it is unknown whether or not they are Green's 
functions. The proof should be performed as to whether 
they satisfy the conditions of Green's functions before 
adapt ion of Green's functions. Eq.(3) can be written in 
component form as 



^ikq^klm 



d 2 E 1 



+ -~Hr{er)ijEj = -iujflQfirJi- (A2) 



By inserting Eq.(Al) into Eq.(A2), one can calculate 
the next equation: 



/ 



X^ikq^klr 



d 2 GE ls + I 

dx doCq C* 



: p r (e r ) tj GE js }U s dv' = U l , 
(A3) 

where Ui = —ico^o^rJi for the source term. Inside the 
bracket is V x V x GE — k^\i T ^ T GE, and to maintain the 
equality, as expected, the equation arrives at the defini- 
tion of the Green's functions: 



V x V x GE - fearer 



d 2 GE 



= IS(r - r') . (A4) 



GE and I are 3-by-3 matrices and the above equation is 
simply a definition of the Green's function. 

The other tensor Green's functions used in this article 
all have similar forms. Moreover, it is easy to prove that 
they meet the definition of the Green's function. The 



G 



equation in Eq.(A4) is essentially nine equations for GEij become decoupled and much simpler to solve through 
that are coupled to each other. However, the equations diagonalization of e r . 
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